Odd-frequency Pairs and Josephson Current through a Strong Ferromagnet 
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We study Josephson current in superconductor / difltusive ferromagnet /superconductor junctions 
by using the recursive Green function method. When the exchange potential in a ferromagnet 
is sufficiently large as compared to the pair potential in a superconductor, an ensemble average of 
Josephson current is much smaller than its mesoscopic fluctuations. The Josephson current vanishes 
when the exchange potential is extremely large so that a ferromagnet is half-metallic. Spin-flip 
scattering at junction interfaces drastically changes the characteristic behavior of Josephson current. 
In addition to spin-singlet Cooper pairs, equal-spin triplet pairs penetrate into a half metal. Such 
equal-spin pairs have an unusual symmetry property called odd-frequency symmetry and carry the 
Josephson current through a half metal. The penetration of odd-frequency pairs into a half metal 
enhances the low energy quasiparticle density of states, which could be detected experimentally by 
scanning tunneling spectroscopy. We will also show that odd-frequency pairs in a half metal cause 
a nonmonotonic temperature dependence of the critical Josephson current. 

PACS numbers: 74.50.-|-r, 74.25.Fy,74.70.Tx 



I. INTRODUCTION 

Ferromagnetism and spin-singlet superconductivity 
are competing orders because the exchange potential 
breaks down spin-singlet pairs. Spin-singlet pairs, how- 
ever, do not always disappear under the influence of 
an exchange potential. Long time ago Fulde-Ferrell^ 
and Larkin-Ovchinniko\*^ discussed inhomogeneous spin- 
singlet superconductivity in the presence of an exchange 
potential. It was shown that the superconducting or- 
der parameter oscillates in real space because the ex- 
change potential shifts the center-of-mass momentum of 
a Cooper pair. Similarly, a Cooper pair has been dis- 
cussed in superconductor / ferromagnet (SF) and super- 
conductor / ferromagnet / superconductor (SFS) junc- 
tions^iii^i^i^i^iiS. Thcsc studies showed that a pairing 
function in a ferromagnet changes its sign periodically 
in real space. As a consequence, SFS junctions may un- 
dergo so-called O-tt transition with varying length of a 
ferromagnet or temperature. 

Previous theoretical studies of the proximity effect in 
a ferromagnet were mainly based on solving the qua- 
siclassical Usadel equations!^ valid when the exchange 
potential Vex is comparable to or smaller than the pair 
potential in a superconductor at zero temperature Aq. 
Cooper pairs can penetrate into a ferromagnet within a 
short distance = \/D/Vex, where D is the diffusion 
constant in a ferromagnet. Thus, penetration of spin- 
singlet Cooper pairs into a ferromagnet with large Vex 
would be impossible and the Josephson coupling via such 
a strong ferromagnet would be vanishingly small. A re- 
cent experiment^, however, demonstrated the existence 



of Josephson coupling through a strong ferromagnet with 
Vex ^ Aq. In addition to this, the experiments^ has even 
shown Josephson coupling in superconductor /half metal 
/ superconductor (S/HM/S) junctions. A half metal is 
an extreme case of a completely spin polarized material 
because its electronic structure is insulating for one spin 
direction and metallic for the other. Thus one has to seek 
a new state of Cooper pairs in a strong ferromagnet. The 
experiment by Keizer et. al. has motivated a number of 
theoretical studies in this directio n -'^'^'^^i^^'^^ . 

Prior to the experiment Eschrig et.al}^ have ad- 
dressed this challenging issue. In the clean limit, they 
have shown that p-wave spin-triplet pairs induced by 
spin-flip scattering at a junction interface can carry 
Josephson current. In practical S/HM/S junctions, how- 
ever, a half metal is close to the dirty limit in the diffu- 
sive transport regime; the elastic mean free path £ may 
be smaller or comparable to the coherence length and is 
much smaller than the size of the half metal L^- Thus, 
the effects of the impurity potential on the Josephson 
current should be clarified in a SFS junction consisting 
of a strong ferromagnet. In this paper, we discuss the 
Josephson effect in SFS junctions for arbitrary magni- 
tude of Vex- When Vex is much larger than Aq, an en- 
semble average of the Josephson current is much smaller 
than its mesoscopic fiuctuation o"'^^i^° . Fluctuations of the 
pairing function in a ferromagnet is responsible for the 
large fluctuations of Josephson current. The Josephson 
current vanishes in S/HM/S in the absence of spin-flip 
scattering at junction interfaces. Spin-flip scattering at 
junction interfaces drastically changes the characteristic 
behavior of Josephson current and properties of Cooper 
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pairs in a ferromagnet. Spin-flip scattering allows for 
the penetration of equal-spin-triplet Cooper pairs which 
have unusual symmetry property called odd-frequency 
symmetry--. When the contribution of equal-spin-triplet 
Cooper pairs to the Joscphson current is dominant, the 
self-averaging property of the Josephson current is re- 
covered. In particular in diffusive S/HM/S junctions, all 
Cooper pairs in a half metal are in the odd-frequency 
equal-spin-triplet pairing statei^. We also discuss local 
density of states in a ferromagnet which reflects the ex- 
istence of odd-frequency Cooper pairs. A part of this 
study has been already published elsewhere^^^. Through- 
out this paper, we use the unit of h — — ^, where ks 
is the Boltzmann constant. 

This paper is organized as follows. In Sec. II, we ex- 
plain the model of SFS junctions on two-dimensional 
tight-binding lattice and the method of calculation. The 
characteristic features of Josephson current in SFS junc- 
tions are discussed in Sec. III. In Sec. IV, we intro- 
duce spin-flip scattering at junction interfaces and dis- 
cuss symmetry properties of Cooper pairs in a ferromag- 
net. We propose an experiment to observe odd-frequency 
pairs in SFS junctions based on calculated results of local 
density of states in Sec. V. The conclusions are formu- 
lated in Sec. VI. 



II. MODEL 



hopping integral t is considered among the nearest neigh- 
bor sites. In a ferromagnet, on-site potential is given 
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FIG. 1: (Color online) (a) A schematic figure of a SFS junc- 
tion on tight-binding lattice, (b) Density of states for each 
spin direction. The Josephson junction is of the SNS, SFS, 
and S/HIVI/S type for Vex/t = 0, 1 and 2.5, respectively. 



Let US consider the two-dimensional tight-binding 
model as shown in Fig. [H^a). A vector r — jx + my 
indicates a lattice site, where x and y are unit vectors in 
the X and y directions, respectively. A junction consists 
of five segments: a ferromagnet (3 < j < Ljv — 2), two 
thin ferromagnetic layers (j = 1, 2, L^ — l, and Ljv), and 
two superconductors (— oo < j < and Lat-I-I < j < oo). 
In the y direction, the number of lattice sites is W and we 
assume a periodic boundary condition. Electronic states 
in a superconducting junction are described by the mean- 
field Hamiltonian 
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-V{r)-&dry, (2) 
Ary =eyHAa2 5ry, (3) 
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(Cj, is the creation (annihilation) operator 
of an electron at r with spin ct = ( t or | ), S in the 
summation means superconductors, aj with j =1 - 3 are 
the Pauli matrices, and (Tq is the 2x2 unit matrix. The 



randomly in a range of — V7/2 < < V//2, where we 
take the probability distribution for e^. unifom on this 
interval, and Cr at different points are uncorrelated. The 
uniform exchange potential in a ferromagnet is given by 
V{r) — Vex^s, where ej for / = 1 — 3 is a unit vector in 
spin space. The Fermi energy /i is set to be 2t in a normal 
metal with Vex = 0, while a ferromagnet and a half metal 
are respectively described by Vex/t = 1 and 2.5 as shown 
in Fig.[Hb). Spin-fiip scattering is introduced at j = 1, 2, 
Ljv — 1, and Ln, where we choose V{r) = Vse2- In a 
superconductor, we take er- = and A is an amplitude of 
the pair potential in s-wave symmetry. The macroscopic 
phases are given by (fj = in the left superconductor 
and by ^pj = ipji in the right one. 

The Hamiltonian is diagonalized by the Bogoliubov 
transformation, 
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where ^\ ^ (7a, ct) is the creation (annihilation) operator 
of a Bogoliubov quasiparticle. The wave functions, u\ 
and wa, satisfy the Bogoliubov-de Gennes cquation^^. 
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The eigen value matrix E\ is diagonal and depends on 
spin channels. To solve the Bogoliubov-de Gennes equa- 
tion, we apply the recursive Green function metho d^^'^'^ . 
In this method, we calculate the Matsubara Green func- 
tion 
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^u;^ + Ex]-'[vi{r'),u{ir')], (8) 
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where a;„ = (2n -I- 1)ttT is a Matsubara frequency, n 
is an integer number, and T is the temperature. The 
Josephson current is given by 

w 



with r' — r + X. In this paper, 2x2 and 4x4 matrices 
are indicated by • " • and ■, respectively. 

In simulations, we first compute the Josephson current 
for a single sample with a specific random impurity con- 
figuration. After calculating the Josephson current over 
a number of different samples, an ensemble average of 
the Josephson current and its fiuctuations are obtained 
as 



6J=^{J^)-{JY 



(11) 
(12) 



where Ji is the Josephson current in the i th sample 
and Ns is the number of samples. Strictly speaking, Ns 
should be taken to be infinity. In this paper, we increase 
Ns until sufficient convergence of (J) and 5J is obtained. 
In the following, Ns is typically taken to be 100-2000. 

To study the characteristics of Cooper pairs in a ferro- 
magnet, we also analyze the anomalous Green function in 
Eq. The pairing function is defined by the anomalous 
Green function and is decomposed into four components. 
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where r = jx + my, lOc = O.OIAq is a low energy cut-off 
and the pairing functions are averaged over whole lattice 
sites at j before ensemble averaging. In Eq. (fT^ . fnifs) 
is the pairing function of spin-singlet (spin-triplet) pairs 
with spin structure of (|ti) - (-1-) \i1))/V2. The pairing 
functions of ItT) and | JJ.) pairs are given by = i/2 — /i 
and fii = i/2 + fi, respectively. 

The quasiclassical Green function methodiii^ is a 
powerful tool to study the proximity effect when the pair 
potential is much smaller than the Fermi energy. The 



quasiclassical Green function, however, cannot be con- 
structed in a half metal because the Fermi energy for one 
spin direction is no longer much larger than the pair po- 
tential. On the other hand, there is no such difficulty in 
our method. These are advantages of the recursive Green 
function method. Throughout this paper we fix the fol- 
lowing parameters: = 74, W = 25, /x = 2t, and 
Vi — 2t. This parameter choice corresponds to the diffu- 
sive transport regime in the N, F and HM layers^^. The 
results presented below are not sensitive to variations of 
these parameters. 



III. SFS JUNCTION WITHOUT SPIN-ACTIVE 
INTERFACE 

In this section, we do not consider spin-flip scattering 
at the interfaces, (i.e., Vs = 0). We first discuss the 
Josephson current in SFS junctions as shown in Fig. [2l 
where T = O.lTc, Aq = O.Oli, Tc is the superconducting 
transition temperature, and the phase difference across a 
junction (p = ipL — is fixed at 7r/2. The presented re- 
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FIG. 2: Josephson current is plotted as a function of the ex- 
change potential Vex for Vs = Q and Ao = O.Qlt. At Vex = 2t, 
a ferromagnet becomes half-metallic as indicated by an arrow. 
The vertical axis is normalized by Jo which is an ensemble av- 
erage of Josephson current at Vex = 0. The number of samples 
used for averaging Ns is 500. 

suits are normalized by Jo which is the ensemble averaged 
of Josephson current in the superconductor / normal 
metal / superconductor (SNS) junctions (i.e.. Vex = 0). 
The Josephson current oscillates as a function of Vex and 
changes its sign almost periodically. The sign changes of 
(J) correspond to the O-tt transition of a SFS junction. 
At the same time, the amplitude of (J) decreases rapidly 
with increasing Vex- For Vex > O.lt, we should pay at- 
tention to the relation (J) <C SJ which means that the 
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Josephson current is not a self-averaging quantity. It is 
impossible to predict the Josephson current in a single 
sample Ji from (J) because Ji strongly depends on the 
microscopic impurity configuration. In fact, the Joseph- 
son current flows in a single sample even if (J) = at 
the transition points. Roughly speaking, (J) vanishes 
because half of samples are 0-junctions and the rest are 
the TT-junctions^*''^^. Since (J) = 0, 5, J approximately 
corresponds to the typical amplitude of the Josephson 
current expected in a single sample. In Fig. [21 we also 
show (I J|), which agrees well with 5 J even quantitatively. 
The relation ( J) = has different meaning for SFS and 
S/HM/S cases. In a SFS junction, the fact that (J) = 
at the transition points is a result of ensemble averag- 
ing and the Josephson current remains finite in a single 
sample. The characteristic temperature and length of a 
ferromagnet at the — tt transitions vary from one sample 
to another. In S/HM/S junctions at Vex = 2.5t, however, 
( J) = means that the Josephson current vanishes even 
in a single sample because (5 J = holds at the same time. 
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FIG. 3; Pairing functions are plotted as a function of position 
j for (a) SNS at Ve:,/t = and (b) S/HM/S at Ve^/t = 2.5, 
where Vs = 0, Ao = O.Olt and Ns = 200. The vertical axis is 
normalized by a pairing function in a superconductor /s. 




FIG. 4: Pairing functions in SFS junctions with Vcx/t — 1, 
Vs = 0, Ao = O.Olt and iV, =500. 



The large fluctuations of Josephson current were dis- 
cussed by Zyuzin et. ali^ by using the diagrammatic 
expansion. An ensemble average of critical Josephson 



current and its fluctuations have a relation for W ^ L 
1 exp(-LAr/fT) 



N 



SJ 

(J) 



(14) 
(15) 



where Exh = D/Lj^ is the Thouless energy, = 
yjDjlixT, and — \/D/Vx is estimated to be about 
four lattice constants (See also Appendix A). In the sec- 
ond line, we replace T by Ag because a measuring tem- 
perature must be smaller than Tj.. For a weak ferromag- 
net (i.e.. Vex ^ Aq), the ratio can be less than unity and 
the Josephson current is a self-averaging quantity. On the 
other hand, in a strong ferromagnet with Vex ^ Ao, the 
ratio becomes much larger than unity. Thus the large 
fluctuation of Josephson current is a robust feature of 
SFS junctions with 14a; 3> Aq. The only way to suppress 
fluctuations is taking the junction width sufficiently large 
because fluctuations are a mesoscopic effect. 

The origin of the large fluctuations in the Josephson 
current can be understood by analyzing pairing functions 
of Cooper pairs. We plot a pairing function of spin-singict 
pairs /o in an SNS junction as a function of position in 
a normal metal j in Fig. [3]^ a), where j = 1 and 37 cor- 
respond respectively to the junction interface and the 
center of the normal metal. The pairing function is cal- 
culated for if = Q and is normalized by its bulk value 
in a superconductor Jb- In SNS junctions, (/o) is a real 
value and is almost constant as shown in Fig.[3I^a), which 
means that spin-singlet Cooper pairs exist everywhere in 
the normal metal. The pairing function for spin-singlet 
pairs in SFS junctions is shown in Fig. [3)Ja). An av- 
erage (/o) decreases exponentially with j according to 
exp(— j/,^;i). At the same time, (/o) oscillates in real 
space and changes its sign. In addition to spin-singlet 
pairs, opposite-spin-triplet pairs appear in a ferromag- 
net for Vex 7^ 0. Since /a is a pure imaginary value, the 
imaginary part of (/a) is plotted in Fig. [DJb). The be- 
havior of (/a) is qualitatively the same as that of (/□) in 
Fig-HUa)- Thus opposite-spin-triplet pairs also contribute 
to the Josephson current. Both 5Jq and SJ^ remain fi- 
nite at the center of a ferromagnet j = 37. Spin-singlet 
and opposite-spin-triplet pairs penetrate deeply into a 
ferromagnet far beyond S^h even though (/o) and (/a) are 
almost zero there. We numerically confirm the relation 
(5/o oc e"-'/'^^, in agreement with Ref. 20. 

In Fig. [5] (a) and (b), we show /o and /a in SFS junc- 
tions for three samples with different impurity distribu- 
tion. The vertical axis is shifted as indicated by hori- 
zontal lines. The pairing functions are in phase near the 
interface [j < C/i)? whereas they are out of phase far from 
the interface. Although the pairing function in a sample 
has a finite value for j > £,h, an ensemble average of them 
vanishes. Cooper pairs do exist in a single sample of fer- 
romagnet even for j ^ ^h- Mesoscopic fluctuations of 
the pairing function provide the origin of the large fluc- 
tuations in the Josephson current. In S/HM/S junctions. 
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as shown in Fig.[3Ub), (/o) and (/s) vanish for j I. We 
have also confirmed that Sfo = Sf^ = for j ^ 1 at the 
same time. Thus, no Cooper pairs exist in a half metal 
for Vs = 0. 
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FIG. 5: Pairing functions in three different samples of SFS 
junction at Vex/t = 1. 

Since (J) <C SJ, the temperature dependence of 
Josephson current also depends on the impurity configu- 
ration. In Fig. [6l we show the Josephson critical current 
as a function of temperature for five different samples, 
where the critical current is estimated from the current- 
phase relation at each temperature. The solid line in 
Fig. [DJa) corresponds to a SFS junction in the 0-state, 
where the critical current monotonically increases with 
the decrease of temperature. On the other hand, the 
broken line corresponds to a junction in the 7r-state. In 
Fig. [6ljb), a junction undergoes the transition from to 
TT state when temperature decreases across 0.5Tc. On the 
contrary, the 0— state is more stable than the tt— state at 
low temperatures in Fig. [H^c). The Josephson current 
is decomposed into a series of J = J2k=i sin(A:(/?). In 
Fig. [SI Ji = characterizes the O-tt transition tempera- 
ture. At the transition temperature, the critical current 
is not exactly zero because a higher harmonic such as 
J2 sin(2(^) contributes to the Josephson current. Some 
SFS junctions undergo the O-tt transition twice as shown 
in Fig.[5fd). The temperature dependence of the critical 
current in one sample can be very different from that in 
another samples. 



IV. SFS JUNCTION WITH SPIN-ACTIVE 
INTERFACE 

The relation (J) <^ 5, J is a. characteristic feature of 
the Josephson current in diffusive SFS junctions with 
Vex ^ Aq. This feature, however, is drastically changed 
by spin-flip scattering at junction interfaces. In this sec- 
tion, we study the Josephson current in the presence of 
spin-flip scattering, (i.e., Vs 7^ 0). In Figs. [7] (a) and (b), 
we show (J) as a function of the spin-flip potential Vs 
for Vex/t = 1 and 2.5, respectively. In both cases (a) and 
(b), we find that |(J)| > 6J for Vs > OM. The Joseph- 
son current recovers the self-averaging property in the 
presence of spin-flip scattering. Reasons can be found 
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FIG. 6: Critical current versus temperatures for five differ- 
ent samples of SFS junction, where Vex/t = 1, Vs = 0, and 
Ao/t = 0.01. 



o.so 



0.00 ijo=,c.^="°0»a0<,aro„<^. 



y -O.SD - 



-1.00 



T,T,=0,1 

I I . I i I I L_ 



0.1P 
A O.OO- 



-0.10 



< • 

I I I . I I ! : 



0.0 0.2 0.4 (».fi Q.S VO OlO 0.2 O.i O.E OjS 1.0 

V.Ji V.J! 



FIG. 7: (a) Josephson current and its fluctuations for T = 
O.lTc, (f = ■k/2 and Ns — 200 as a function of interface spin- 
flip scattering Vs for Vex/t = 1 and (b) for Vexjt = 2.5. 
The vertical axis is normalized by an ensemble average of 
Josephson current at Vex ~ and Vs = 0. 



by analyzing the pairing functions in a ferromagnet, as 
shown in Figs. [S] and [51 where four pairing functions are 
plotted as a function of position j at Vs/t = 0.4. In SFS 
junctions as shown in Fig.[5]Ja), equal-spin-triplet Cooper 
pairs penetrate into a ferromagnet by spin-flip scattering 
at interfaces. Although averages of the pairing function 
for opposite-spin pairs vanish at j ~ 37, their fluctua- 
tions remain finite as shown in Figs.[5Ja) and (b). Thus 
four types of Cooper pairs carry the Josephson current 
in a SFS junction. In a S/HM/S junction, on the other 
hand, only "ff-pairs exist in a half metal as shown in 
Fig. [Hfa) and (b). The pairing functions (fn), (/o), and 
if 3) vanish for j 1. We note that fluctuations of these 
pairing functions behave similar to their averages. In 
both SFS and S/HM/S, (/|-f) becomes much larger than 
dfo because the exchange potential does not break down 
equal-spin-triplet Cooper pairs and f-\-\ does not suffer 
sign change in real space. Thus the Josephson current 
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becomes a self-averaging quantity as shown in Figs. [Tlja) 
and (b). 
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FIG. 8: Pairing functions in SFS junctions are plotted as a 
junction of j. Ensemble averages and some of their fluctua- 
tions are shown in (a) and (b), respectively. The number of 
samples are taken to be 500. 
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FIG. 9: Pairing functions in S/HM/S junctions are plotted 
as a junction of j. The number of samples are taken to be 
200. 

Here we address an unusual symmetry property of 
Cooper pairs in SFS junctions. In Fig. llOi we show four 
pairing functions in a SFS junction as a function of cj„, 
where j = 37, Vs = 0.2t, = 0, and Vex = t. Al- 
though the Green function at cj^ = is not defined, we 
put /ff = = = at aj„ = 0, and connect results 
for positive ujn with those for negative w„. The pairing 
function /o is an even function of ujn, whereas /||, fn, 
and /a are an odd function of ujn^. Since electrons obey 
Fermi statistics, pairing functions must be antisymmetric 
under interchanging two electrons. 



(16) 



where [/]* denotes the transpose of / meaning the in- 
terchange of spins, ft is well known that ordinary even- 
frequency pairs are classified into two symmetry classes: 
spin-singlet even-parity symmetry and spin-triplet odd- 
parity one. In the former case, the negative sign on the 
right hand side of Eq. p6)) arises due to the interchange 
of spins, while in the latter case due to r ^ r'. In the 
present calculation, all components on the right hand side 
of Eq. (jl3p have s-wave symmetry. This is because the 



pairing functions are isotropic in both real and momen- 
tum spaces due to diffusive impurity scatterings^^. As a 
result, /if, fii, and /a must be an odd function of uJn to 
obey Eq. The fraction of odd-frequency pairs de- 

pends on parameters such as the exchange potential and 
the spin-flip potential. As shown Fig. [3]Ja), all Cooper 
pairs have even-frequency symmetry in SNS junctions at 
Vex = and Vs = 0. Even- and odd- frequency pairs 
have almost same fraction in SFS junctions at Vex — t 
and Vs = as shown in Fig. [H In the presence of spin- 
flip potential, odd-frequency pairs become dominant as 
shown in Fig. [5] In particular, all Cooper pairs have 
odd-frequency symmetry in S/HM/S junctions as shown 
in Fig. [ni The Josephson current in Fig. ^h) is carried 
purely by odd-frequency pairs in S/HM/S junctions. 
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FIG. 10: Pairing functions in a SFS junction are plotted as 
a junction of u!„ for V^x/t = 1 and Vs/t = 0.2. 

The results in Fig. [7] show that the amplitude of 
Josephson current first increases with the increase of Vs 
then decreases. Here we discuss the analytical expression 
of the Josephson current in S/HM/S junction at T = 0, 

(J) = -Ji 

Vl-Vr- V^'^ ^i'^ ) sin + eg • ( Vl X ) cos ^ 

(17) 

(18) 



Ji = ^^eETh gNb^ > 0, 



TT 

1 /■'^/2 



COS''' 7 



'V^ + i) cos* 7 - V^ cos2 -f + Vf 



(19) 



Here Vr = ^Li ^4''^efe and Vl = EL 
the dimensionless magnetic moments at the right and 
left junction interface, respectively. We assume that 
\Vr\ = l^il = Vs and {e^/h)gM is the normal con- 
ductance of a half metal. Details of derivation are dis- 
cussed in Appendices A and B. To compare Ea. p?|) with 
the results in Fig. [Hb), we choose V r — V l = ^562. 
We also note that the magnetic moment in a half metal 
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Vexes- For Vs <C 1, the amplitude of the 



Josephson current increases with Vg because b — 2 and 



Vl ■ Vr — V^^^vj^^ ~ ^s- this case, spin-flip scat- 
tering assists the Josephson current. For large T/g, on 
the other hand, the Josephson current decreases propor- 
tionally to V^^ because the spin-flip potential acts like a 
potential barrier and suppresses the transmission prob- 
ability of the interface. The Josephson current shows 
reentrant behavior as shown in Fig. [71 The Josephson 
current in Fig. [TUb) is calculated at = tt/2. The re- 
sults indicate that the S/HM/S junction is a 7r-junction. 
This conclusion depends on the direction of the mag- 
netic moments at the spin-flip interfaces. In the case of 
Vl = V R, the Josephson current in Eq. (|17p is propor- 
tional to — Ji siniy9 in agreement with Fig. Efb). In the 
case of antiferromagnetic alignment, Vl — —Vr, the 
junction is in the 0-state. Thus we conclude that the 
stability of the 0-state and that of the 7r-state depend 
on the alignment of the magnetic moments at the two 
interfaces22i2^. This feature indicates a new direction 
to controlling the 7r-phase shift by using ferromagnetic 
materials. For Vex ■ {Vl ^ Vr) ^ 0, the Josephson 
current flows even at tp — because such spin configura- 
tion breaks the chiral symmetry of a junction. We have 
numerically confirmed the Eq. (|17[) . 

In the end of this section, we discuss the temper- 
ature dependence of the Josephson critical current in 
SFS and S/HM/S junctions in Fig. [Til where we choose 
Ao = O.OOSi in connection with the density of states in 
the next section. The Josephson current has almost a si- 
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FIG. 11: Temperature dependence of critical Josephson cur- 
rent in SFS (a) and S/HM/S (b) junctions. Ns is taken to be 
200. 

nusoidal current-phase relationship. The critical current 
for Vs/t — 0.2 and 0.3 in a SFS junction first increases 
with the decrease of temperature then decreases as shown 
in Fig. Illf a). Such reentrant behavior is seen more 
clearly in a S/HM/S junction as shown for Vs/t — 0.2, 
0.3 and 0.4 in Fig. fTlT b). In a Josephson junction consist- 
ing of conventional s-wave spin-singlet superconductors, 
such reentrant behavior is very unusual. This behav- 
ior has also been reported in Ref. [l^. The results for 
Vs/t = 0.4 in SFS and Vs/t ^ 0.5 in S/HM/S junc- 
tions, on the other hand, show a monotonic temperature 
dependence. 



V. DENSITY OF STATES 

The proximity effect changes the low energy spectra of 
a quasiparticle in a normal metal. In a SNS junction, it is 
well known that the penetration of usual even-frequency 
spin-singlet s-wave Cooper pairs suppresses the quasipar- 
ticle density of states for E < Eth- This suppressed den- 
sity of states is called minigap. In this section, we discuss 
the proximity effect of odd-frequency pairs on the quasi- 
particle density of states. In our method, the density of 
states is given by 



N{E,j) = -iip ImTrG£+,^(r,r), (20) 



w 



ttW 



m—l 



where 7 is a small imaginary part. In Fig. 1121 we show 
the local density of states (LDOS) at j = 37 in SFS 
junctions, where ip — 0, and 7 = O.IAq. The results for 
S/HM/S junctions are presented in Fig. [HI The LDOS 
is normalized by its value at E = 1.2Ao. Here we choose 
Ao = O.OOSt so that Eth ^ 0.3Ao is slightly smaller than 
Aq. In the absence of spin-flip scattering, the ensemble 
average of LDOS is almost constant in both Figs. [T^ and 
[T3I At Vs/t = 0.3, the penetration of odd-frequency 
pairs enhances LDOS for E < O.SAq. On the other 
hand, LDOS is suppressed around E ^ O.8A0 because of 
a sum rule for the density of states (i.e., / dE N{E,j) = 
const.). The low energy spectra of LDOS increase with 
increasing Vs/t as shown in Figs. [T^ b) and [T^Tb). At 
Vs/t = 0.5, LDOS has a peak at E = 0. Thus the pene- 
tration of odd-frequency pairs enhances the quasiparticle 
density of states for E < Etii- This tendency is just 
opposite to the minigap structure due to penetration of 
even- frequency pairs. The shape of the zero-energy peak 
in Figs.fTWb) andfTSTb) is almost independent of the po- 
sition in a half metal. The peak is much stronger than 
the enhancement of the LDOS found in weak ferromag- 
nets^iiSi^i^i^. In such SF junctions^i, the LDOS has an 
oscillatory peak/dip structure at E = 0, which rapidly 
decays with the distance from the SF interface. There- 
fore, the large peak at = in the LDOS is a robust 
and direct evidence of the odd-frequency pairing in a fer- 
romagnet. Scanning tunneling spectroscopy (STS) could 
be used to detect such a peculiar pairing state. 

As shown in Fig. [T21 however, the penetration of odd- 
frequency pairs does not always give rise to a zero-energy 
peak in LDOS. The results for Vs/t = 0.3 and 0.4 have 
a broad peak at a finite energy smaller than Eth. This 
situation is slightly different from the large zero-energy 
peak in a normal metal due to the penetration of odd- 
frequency pairs from spin-triplet odd-parity supercon- 
ductors2^i^i2^i^. In a spin-triplet superconductor junc- 
tion, LDOS in a normal metal always has a large zero- 
energy peak because a midgap Andreev resonant stated 
assists the zero-bias peak. In ferromagnetic junctions, 
on the other hand, such a quasiparticle state is absent 
at the junction interface. The broad peak structure in 
the LDOS is responsible for the nonmonotonic temper- 
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FIG. 12: Local density of states at j = 37 in a ferromagnet 
of SFS junction, where Ns = 2000. 
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FIG. 13: Local density of states at j = 37 in a half metal of 
S/HM/S junction, where N, = 2000. 



ature dependence of the critical current in Fig. [TIT b). 
At high temperatures, quasiparticle states around the 
peak contribute to the Josephson current. At low tem- 
peratures, however, such quasiparticle states cannot con- 
tribute to the Josephson current'^^. We conclude that 
odd-frequency pairs could also be confirmed by measur- 
ing the dependence of the critical current on temperature. 



VI. CONCLUSION 

In conclusion, we have studied the Josephson effect in 
superconductor / diffusive ferromagnet / superconductor 
(SFS) junctions by using the recursive Green function 
method. When the exchange potential in a ferromag- 
net is much larger than the pair potential in a super- 
conductor, the Josephson current is not a self-averaging 
quantity. This is because spin-singlet Cooper pairs pene- 
trating into a ferromagnet far beyond £,h cause large fluc- 
tuations of the pairing function. As a consequence, the 
temperature dependence of the critical Josephson cur- 
rent in one sample can be very different from that in 
another sample. When a ferromagnet is half-metallic, 
the Josephson current vanishes in the absence of spin- 
flip scattering at junction interfaces. Spin-flip scatter- 
ing at interfaces allows equal-spin-triplet odd-frequency 
Cooper pairs to penetrate into a ferromagnet. The ratio 
of odd-frequency pairs to even-frequency ones depends on 
the exchange potential in a ferromagnet and the spin-flip 
potential at interfaces. The Josephson current recovers 



the self-averaging property when the fraction of equal- 
spin-triplet pairs becomes large. In half-metallic SFS 
junctions, all Cooper pairs have odd- frequency symme- 
try. The penetration of odd-frequency pairs enhances low 
energy quasiparticle density of states in a ferromagnet. 
Such low energy spectra could be probed by scanning 
tunneling spectroscopy and determining a nonmonotonic 
temperature dependence of the critical Josephson cur- 
rent. We also discuss a way to realize a 7r-j unction by 
controlling magnetic moments in ferromagnetic layers. 
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APPENDIX A: FLUCTUATIONS OF 
JOSEPHSON CURRENT 

The purpose of this appendix is to explain Eq. (|14p . 
Since fluctuations of Josephson current have been calcu- 
lated by the diagrammatic expansioni2i2£ii2S, we also cal- 
culate the Josephson current in SFS junctions in the same 
method. We assume that relations Exh ^ Aq, Vex ^ M 
and Ln ^ £ are satisfied. In the lowest coupling, the 
Josephson current is given by a formula^ 



J : 



Tr r 



eh Th 



l.r 



(Al) 

where I (r) denotes a propagating channel at the left 
(right) junction interface. In Fig. [Mja), a propagation 
process of the first term in Eq. (|Aip is schematically il- 
lustrated. We calculate the transmission coefficients in 
a ferromagnet such as and t^^ and Andreev reflection 
coefficients at interfaces such as f^^ and ff^ by parts. 
The Andreev reflection coefficients are calculated at an 
ideal NS interface as shown in the left figure of Fig.fHTb). 
The results are given by 



(72- 



^0 



(A2) 
(A3) 



where fi„ = y/u;^ + Aq. The effect of the exchange po- 
tential is considered through transmission coefficients of 
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an electron in a ferromagnet 



equation with appropriate boundary conditions'^ 



''Ir — 



tui) 



(A4) 



The transmission coefficients of a hole are defined in the 
same way by e — *■ /i in the equation above. The trans- 
mission coefficients are represented by the Green function 
as 

t^,((T) =^e^fc<— JJ dyLdyRY;iyn)YiiyL) 

>^Gl^{xR,yR-XL,yL) (A5) 

tl{a) = - »e-*--+^'='--H«,, j j dyLdyRYC{yL)Yr{yR) 

'x G-ur,(^ii'yR'^L,yL), (A6) 

where Yi{y) is a wave function in the y direction and 
vi = ki/m with ki being a wave number in the x direc- 
tion on the Fermi surface in the I th propagating channel. 
In above expression, we have assumed that two ideal lead 
wires are attached to the both sides of a diffusive ferro- 
magnet, and a:^ < and {xr, > Lj^) are taken to be in 
the lead wires. The Green function is given by 



1 



dk e*'= (''-'-') 



(27r)2y VexS+ ^sgn(tj„) 



, (A7) 



where t is the elastic mean free time, = /ni — and 
s — 1(— 1) for a =1 (J,). An ensemble average of trans- 
mission coefficients is calculated by the diagrammatic ex- 
pansion 



Y,{tl{a)tUcj')) 



Ir 



V 



2 /-W r'W 



f I dyLl dyRP^'^' (L + S, yR- ^6, yv, 2t^„), (A8) 



P^- (r,r') = (G-Jr,r')G^'„„(r,r')), 



(A9) 



where P^'^ is the Cooperon propagator which satisfies 
the equation 

= 2nNoS{r - r'). (AlO) 

Since ^ ^ Vex, the diffusion constant D, the Fermi veloc- 
ity vp and the density of states at the Fermi energy A^o do 
not depend on spin directions. In Fig. ll4f c). we illustrate 
the Cooperon and diffuson propagator, where uji — 2t:IT 
is a boson Matsubara frequency. In Fig. [HT d) we show 
two diagrams which contribute to the Josephson current. 
The left (right) diagram in Fig. [TW d) corresponds to the 
first (second) term of Eq. (jAip . Only P^^ and P^'^ con- 
tribute to the Josephson current because the Andreev 
reflection coefficients are off-diagonal in spin space. To 
calculate the Cooperon propagator, we solve the diffusion 



DV^hxir) = Xhx{r), 
^A(r)Uo,L«=0, 
dhxir) 



dy 



0. 



(All) 
(A12) 

(A13) 



y=0,W 



The Cooperon propagator is represented by using wave 
functions and their eigen values of the diffusion equation. 
The results are 



P^-' (r, r'; 2c^„) = 2^iVo E E H^- 

n=lm=0 ^ 

sin(pnx) sin(p„a;0 cos(i^„ty) cos{vrny') 

2{|u;„| - iVexs{l - <5.,.')} + D{pl + vl) 
WK mil 

l/W form = 
2IW form^O, ' 



Vn — J 1 



, (A14) 
(A15) 
(A16) 



By substituting the above results into Eq. (lASp . we arrive 
at 



^(^F,(a)^,^(a))=5iv 

Ir 

ON ^ttNoDW/Ln, 



\sinh?7°'°' 



2\LO„\-2isVc, 
2|^„ 



(T = (T 



(A17) 
(A18) 

(A19) 



where (2e^/h)gN is the conductance of a ferromagnet, a 
denotes the opposite spin state of cr, and integration in 
the y direction in Eq. (|A8p is carried out at S = i/\/2. 
The Josephson current becomes 



( J) =2e 5Ar sin ((5 T ^ ^ 



By substituting equations 



sinh r]° 



Ao 



5Z si 



sinh 7]° 



= 4 



2Ve, 



Ao 



nr. 



e '-h sin , 

ETh V 6i 

^ 2Ao 



(A20) 

(A21) 
(A22) 



into above expression, the Josephson current at T = 
results in 



16V2 AkT -i^ . (Ln 



X sin If. 



Eth 



(A23) 



This expression is also valid for V^x ^ Aq because the re- 
lation Vex 3> Aq was not explicitly used in the derivation. 
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(b) 



> 




> 






r' I- 




r' 


— > — 




■* 





Thus the ratio is described by 




(A25) 



Since P"" ^ ^-Ln/^t ^ temperature dependence of fluc- 
tuations is also expected to be e~^"/'^^. Thus we arrive 
at Eq. (fH)) . In a recent paper, mesoscopic fluctuations 
of the Josephson current were calculated within the qua- 
siclassical Green function technique'*^. In this approach 
the fluctuations are slightly larger than those within the 
diagrammatic expansio n^^i'^^ . The difference may stem 
from the proximity effect on electronic structure in a 
normal metal such as the minigap in the quasiparticle 
density of states. In the diagrammatic expansion, such 
effect is not taken into account. 



p; 



8,1 (J 



(d) 



CD 




FIG. 14: (a): a propagation process of a quasiparticle in a 
SFS junction. We describe the Andreev reflection coeflicients 
in Eq. (IA1[I by those at an ideal normal- metal/superconductor 
interface as shown in (b). (c): Cooperon and difi'uson prop- 
agator. The diagrams for Josephson current (d) and its fluc- 
tuations (e). 



In Fig \Wl e) , we show two typical diagrams for fluc- 
tuations. Not only P°'°' but also P°'°' contributes to 
fluctuations. The Cooperon P'^'^ behaves like e"^"/'''' 
similar to the Josephson current. On the other hand, 
P'^'^ ^ e-L«/iT_ Thus amplitude of fluctuations in SFS 
junctions is almost the same as that in SNS junctions. 
Our approach, however, is not suitable for calculating 
fluctuations because a number of diagrams contributes 
to fluctuations in addition to Fig.fTHe). Here we present 
the result for SNS junctions at T = and Lp ^ tt/2 ob- 
tained by a slightly different approach^ 



6J 



eExh 



(A24) 



The fluctuations in SFS junctions are given by SJ/\/2 
because contribution of P"^'"^ is negligible for Vex ^ Eth- 



APPENDIX B: NEGATIVE JOSEPHSON 
COUPLING 

Here we express the Josephson current in S/HM/S 
junctions with spin-active interface on the basis of the 
diagrammatic expansion. In a half metal, we assume 
that the magnetic moment is parallel to 63. Thus trans- 
mission coefflcients of an electron become 



(Bl) 



because electric structure for [ spin is insulating in a half 
metal. Transmission coeflicients of a hole are defined in 
the same way by e — » ft-. Andreev reflection coefficients 
f^" and ff in Eq. are calculated at a normal- metal/ 
superconductor interface at which a spin-fiip potential 
vpV R-o'5{x— L]si) is introduced as shown in the right fig- 
ure of Fig. [Tif b) . Andreev refiection coefficients f^'^ and 
ff'* are also calculated at a normal- metal/ superconduc- 
tor interface at which a spin- flip potential vp'V l'^5{x) is 

considered. We assume that Vi(fl) = ^L(R)^k and 

\V l\ = \y r\ = Vs- The calculated results of Andreev 
reflection coefficients are given by 



r'r=-Qi^2 [Aido + iBiV L ■ or] 
rt''=Qi [Ai&o + iBiV L ■ a-] aae 



e 



Ao 



1i 



Ai = - n„v^ + qt 

Bl ^qi{fln + UJ„), 



(0„ -I- UJn) 



(B2) 
(B3) 

(B4) 

(B5) 
(B6) 



where qi = ki/kp > are normalized wave number of 
the I th propagating channel in the current direction. 
Andreev refiection coefficients at the right interface f'^'^ 
ff'- are also obtained hy I ^ r, Vl ^ Vr, and 
-!■ ifR in above expression. Since the half metal is in 



and 



the diffusive transport regime, transmission coefficients 
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across the half metal, namely tjj,, and tj?; are in- 

dependent of propagating channels I and r. Thus average 
of the Andreev reflection coefficients over all propagat- 
ing channels contribute to the Josephson current^. We 
define such Andreev reflection coefficients as 



'he 



(B7) 



l{r) 



- CT2 [a&a + ibVuR) ■ a-] g-^'^MH) , (B8) 



l{r) 



(BIO) 



where Nc — Wkp /i^ is the number of propagating chan- 
nels at Fermi energy, a and b are real numbers depending 
only on ujn , Aq , and Vs ■ A part of Eq. (jAip becomes 



r 



Lr 



X {^0 - ^3)ia + ibVL ■ (t)(o-o + 0-3)] , 



(Bll) 



(B12) 



sinhj^TT 

X Vl-Vr~ l^i'V^'^ + ieg • {Vr X Vl)] , (B13) 
where we used Eq. (jA17p . In the same way, we obtain 
l2=J2{Tiiff-i';rf^''-il]), (B14) 



Lr 



9n- 



sinh?7TT 

Vl-Vr- V^i'V^'^ - ie^ ■ {Vr X Vl) 



(B15) 



As a result, the expression for the Josephson takes the 
form 



(J)«-Ji[( 

+63 • {Vl X Vr) cos If] , 



sinhryTT 



b^ > 0. 



(B16) 
(B17) 



The Josephson current is zero in the absence of spin-flip 
scattering at the interface (i.e., Vl = Vr = 0). We note 
that the ratio rj^^ /smb. rj^^ rapidly decreases to zero for 
ujn/ETh ^ 1, whereas Un dependence of b is scaled by 
Aq. For Erh < Ao, we find at T = 



Ji = eErh gNb , 



b = 



7r/2 

c?7 



cos^ 7 



(B18) 



(B19) 



/o (T/|-fi)cos4 7-y|cos2 7 + y4 

Although Eq. (|B16p describes well the dependence of the 
Josephson current on Vl and Vr, it does not explain 
the nonmonotonic temperature dependence of the crit- 
ical current shown in Fig. [TlT b). This is because the 
proximity effect on the density of states in a half metal 
is not taken into account in the above estimate. 
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